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CHAPTER  I 
INTRODUCTION 

In  recent  years  the  investigation  of  the  interaction  of 
fast  heavy  ions  passing  through  matter  has  received  increased 
attention.  The  first  studies  of  electron  charge  exchange 
processes  began  as  early  as  1922,  but  not  until  194$,  when 
accelerators  were  developed,  were  intensified  studies  begun. 
In  particular  the  improved  techniques  for  the  acceleration 
of  heavy  ions,  which  are  still  being  developed  further,  allow 
more  specific  experiments  where  the  nuclear  and  ionic  charge 
of  the  incident  ions  as  well  as  the  energy  can  be  chosen  and 
varied  within  wide  limits.  As  not  too  many  experiments  with 
highly  charged  ions  have  been  done  in  the  past,  this  special 
field  of  heavy  ion  physics  is  potentially  rich  and  may  con- 
tribute not  only  to  an  improvement  of  future  heavy  ion  ac- 
celerators, which  are  also  of  great  use  in  many  other  fields, 
but  in  addition  to  a  better  understanding  of  the  complicated 
charge  exchange  collisions,  a  many  body  problem  which  is 
difficult  to  describe  by  theoretical  methods. 

In  this  work,  a  few  contributions  of  an  auxiliary  but 
basic  character  are  presented.   In  part  II,  the  process  of 
electron  loss  and  capture  is  described,  the  basic  differential 
equations  for  the  charge  exchange  collisions  are  discussed 
and  the  experimental  and  theoretical  situation  prior  to  1969 
is  briefly  summarized.   Parts  III  and  IV  give  two  basic  com- 
puter programs  which  are  usually  of  great  importance  when 


experimental  data  of  heavy  ion  charge  states  are  to  be 
analyzed;  in  the  program  NEDT,  methods  for  the  numerical 
integration  of  the  fundamental  system  of  differential  equa- 
tions for  ionic  charge  states  are  given,  and  in  the  routine 
EQUDIS,the  equilibrium  charge  fractions  are  calculated.   Both 
programs  can  be  used  for  any  sets  of  given  ionic  charge  ex- 
change cross  sections.   In  order  to  use  those  routines  for 
iteration  procedures,  special  attention  was  given  to  keeping 
the  computational  execution  time  at  a  minimum.   Finally,  in 
part  V,  an  analysis  of  the  shape  of  equilibrium  distributions 
is  attempted.   Experimental  distributions,  measured  at  Oak 
Ridge,  have  been  fitted  to  Gaussian  and  more  general  distri- 
butions. Characteristic  deviations  from  a  Gaussian  distribu- 
tion as  well  as  asymmetries  have  been  found.   Distributions 
of  Br  ions  at  20  MeV  and  I  ions  at  15  and  60  MeV  stripped  in 
dilute  gases  of  H2,  He,  N2,  He,  Ar  and  Kr  have  been  analyzed 
in  this  work. 


CHAPTER  II 
DESCRIPTION  OF  THE  CHARGE  EXCHANGE  PROCESS 

When  ions  enter  into  a  gaseous  or  solid  target,  col- 
lisions take  place  between  the  ions  and  the  target  atoms. 
In  such  a  collision,  the  incident  ion  may  either  capture  or 
lose  one  or  more  electrons.   Therefore,  one  can  define  a 
cross  section  <?"ij,  which  describes  the  probability  that  the 
charge  i  of  an  ion  changes  into  j  during  an  encounter.  The 
cross  sections  are  usually  given  in  cm2/atom,  but  sometimes 
the  units  cm2/molecule  are  preferred,  especially  if  the  tar- 
get consists  of  molecules. 

We  consider  only  those  collisions  in  which  the  charge  of 
the  ions  is  actually  changed.  The  cross  sections  introduced 
above  do  not  take  into  account  elastic  collisions  or  encounters 
in  which  only  excitation  occurs.  If  we  assume  further,  that 
the  total  energy  loss  of  all  ions  in  the  target  is  negligible, 
the  relative  fractions,  Y±   ,  of  ions  carrying  a  charge  i  in 
a  beam,  change  with  increased  target  thickness  x,  measured 
in  atoms  of  molecules  per  cm2  as  follows: 


In  this  system  of  equations  all  possible  charge  transfers  are 
taken  into  account.  If  we  neglect  the  probabilities  for  mul- 
tiple electron  exchange  and  consider  only  collisions  in  which 


the  charge  of  the  ions  charges  by  one  unit  of  charge 
(  fj  -  ij  »  1),  Eq.  (1)  reduces  to 

gi  -  ^(i-ur,.,  -  tfyi)  +<rL(i)ii  +n<1+1)Vi  ■    (2.2) 

where  Q£l*>  £l-i   and  (^CD^C-,;^     are  used  for  a  better 
understanding  of  the  influence  of  capture  and  loss  processes 
on  the  charge  state  populations,  as  expressed  in  Eq.  (2). 

The  reduced  form  in  Eq.  (2)  may  be  used. in  connection 
with  light  ions  stripped  in  light  targets  such  as  H2  or  He. 
In  this  case,  the  probabilities  for  multiple  electron  exchange 
are  very  small  and  can  be  neglected.  But  if  the  ions  have  a 
large  atomic  number,  or  if  the  target  is  heavy,  simultaneous 
capture  or  loss  of  more  than  one  electron  in  a  single  col- 
lision has  an  increased  probability  and  cannot  be  neglected. 
Thus  the  system  (1)  is  to  be  used.  For  example,  if  Br  ions 
pass  through  an  argon  target,  double  loss  is  approximately 
half  as  likely  as  single  loss.  As  the  present  work  is  mainly 
concerned  with  stripping  of  heavy  ions,  the  complete  system 
(1)  will  be  considered  in  this  and  the  following  chapters. 

The  range  of  different  charge  states  i  in  Eq.  (1)  should 
include  all  significant  charge  states.   In  principle,  all 
charges  ranging  from  i  =  -1  to  i  =  Z,  the  nuclear  charge  of 
the  ion  in  question,  are  possible.   In  practice,  however,  an 
ion  beam  shows  only  a  limited  number  of  charges  even  after 


having  passed  through  a  thick  target.   For  Br  and  I  ions, 
usually  10  to  12  adjacent  charge  fractions  can  be  observed. 

If  an  incident  beam  consists  of  ions  nhich  carry  the 
same  charge  i,  a  whole  distribution  of  charge  states  will  be 
produced  according  to  the  various  charge  changing  probabili- 
ties Q±y     Thus  the  various  charge  fractions  Y..  will  change 
with  x,  the  thickness  of  the  stripper.   The  corresponding 
distribution  is  called  a  no n- equilibrium  distribution,  which 
is,  of  course,  a  solution  of  Eq.  (1).  An  example  for  such  a 
distribution  is  shown  in  Plates  la  through  c  and  Table  la 
through  c.  The  distribution  depends  on  the  initial  charge 
state  of  the  incident  ions.  As  the  thickness  of  the  stripper 
is  increased  more  and  more,  all  charge  fractions  approach  a 
constant  equilibrium  value.   These  equilibrium  fractions  are 
no  longer  dependent  on  the  special  initial  conditions  or  on 
the  target  thickness,  as  soon  as  the  equilibrium  thickness  is 
reached.   Examples  for  equilibrium  distributions  are  discussed 
in  part  V.  The  equilibrium  distribution  is  determined  by  the 
species  of  ions,  their  energy,  and  by  the  target  material. 
Experiments  have  shown  that  the  energy  loss  is  indeed 
very  small,  as  long  as  x  does  not  greatly  exceed  the  equilib- 
rium value.  Therefore,  one  can  assume  that  there  is  no 
energy  dependence  of  the  cross  sections  within  a  given  system 
(1),  and  one. then  has  constant  coefficients.  Another  effect, 
however,  can  cause  serious  complications.  If  the  target  is 
very  dense,  excited  ions  will  not  have  fallen  back  into  the 
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Table  1-a  Non-equilibrium  charge  distribution  obtained  with 
program  NEDT.   Initial  charge  state  is  1+. 


Ionic  Charge 

Distribution       Non-equilibrium      Initial  Charge  1+ 

X     ZQX    Y(l)    Y{2)    Y(3)   Y(4)     Y(5)    Y(6) 


0.0010  1.00  98.900  1.100  0.000  0.000  0.000  0.000 

0.0020  1.02  97.319  2.165  0.016  0.000  0.000  0.000 

0.0030  1.03  96.752  3-206  0.042  0.000  0.000  0.000 

0.0040  1.04  95.696  4-224  0.078  0.001  0.000  0.000 

0.0050  1.05  94-654  5.220  0.124  0.002  0.000  0.000 

0.0060  1.07  93.623  6.194  0.179  0.003  0.000  0.000 

0.0070  1.03  92.605  7.146  0.243  0.005  0.000  0.000 

0.0080  1.09  91.559  8.073  0.316  0.007"  •  0.000  0.000 

0.0090  1.10  90.604  8.988  0.397  0.010  0.000  0.000 

0.0100  1.11  89.621  9.379  0.485  0.014  0.000  0.000 

0.0200  1.22  80.393  17-751  1-750  0.102  0.004  0.000 

0.0300  1.32  72.192  23.923  3-553  0.312  0.018  0.001 

0.0400  1.42  64.898  28.633  5-695  O.67O  0.052  0.003 

0.0500  1.52  58.402  32.270  8.021  1.184  0.115  0.008 

0.0600  1.62  52.611  34-889  10.417  I.S48  °-217  °-019 

0.0700  1.72  47-443  36.710  12.795       2.652  0.363  0.0J7 

0.0800  1.81  42.824  37-878  15-095  3-577  0.5o0  0.066 

0.0900  1.91  38.689  33.515  17.272  4- 605  0.812  0.108 

0.1000  2.00  34-976  38.726  19.300       5.714  1-113  0.166 

0.2000  2.31  13.417  29.926  29.902  17.670  6.903  2.132 

0.3000  3-46  5. 761  19.021  28.221  24-872  14-719  7-406 

0.4000  3.95  2.751  11-32,1  22.973  26.643  21.138  14-605 

0.5000  4.32  1.450  7-579  17.926  25.573  25.535  21.936 

0.6000  4-53  0.836  5-077  14.018  23.594  28.203  28.273 

0.7000  4.76  0.524  3.583  11.223  21.636  29.736  33.248 

0.3000  ;,.39  0.353  2.670  9.291  20.027  30.724  36.934 

0.9000  4.97  0.255  2.098  7-977  18.812  31-289  39-563 

1.0000  5.03  0.196  1.732  7-036  17.934  31.637  41.414 

2.0000  5.15  0.099  1-082  5-391  16.135  32.220  45-073 

3.0000  5.15  0.098  1.074  5-371  16. ILL  32.226  45-117 

Equilibrium  Distribution 

3.3125  5.15  0.098  1.074  5.371  16.113  32.227  45.117 

Exact  Equilibrium  Distribution  as  Follows 

5-15  0.098  1.074  5.371  16.113  32.227  45-117 
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Table  1-b  Non- equilibrium,  charge  distributions  obtained  with 
program  NEDT.   Initial  charge  state  is  4+- 

Ionic  Charge 

Distribution       Non- equilibrium     Initial  Charge  4* 

X     ZQX    T(l)    Y(2)    Y(3)    T(4)    Y(5)    Y(6) 


0.0010  4.00  0.000  0.000  0.300  93.900  0.800  0.000 

0.0020  4.01  0.000  0.001  0.590  97.226  1.574  0.003 

0.0030  4.01  0.000  0.002  0.874  96.770  2.331  0.022 

0.0040  4-02  0.000  0.004  1.152  95-731  3-072  0.040 

0.0050  4-02  0.000  0.007  1-424  94.7C9  3.797  0.064 

0.0060  4.O3  0.000  0.010  1.689  93-703  4-506  0.093 

0.0070  4-03  0.000  0.013  1.949  92.713  5-199  0.126 

0.0030  4.04  0.000  0.017  2.204  91-733  5.377  0.164 

0.0090  4.04  0.000  0.022  2.452  90.779  6.541  0.206 

0.0100  4.05  0.000  0.027  2.696  39.335  7-190  0.253 

0.0200  4.10  0.001  0.096  4.352  31.170  12.951  0.031 

0.0300  4.14  0.002  0.195  6.563  73-763  17-544  1-929 

0.0400  L.19  o.oozt  0.312  7.911  67.423  21.191  3-159 

0.0500  4.23  0.007  0.440  3.964  6I.963  24.073  4.552 

0.0600  4.27  0.011  0.573  9.774  57.249  2o.340  6.053 

0.0700  4.31  0.016  0.707  10.333  53-162  23.110  7-613 

0.0300  4-35  0.022  0.336  10.343  49.603  29-432  9.214 

0.0900  4-39  0.023  0.961  11.163  46.491  30.536  10.816 

0.1000  4.42  0.035  1.073  11.339  43.756  31.333  12.404 

0.2000  4.70  0.107  1.792  10.353  23.566  33-007  25-669 

0.3000  4-36  0.150  1.392  9-262  22.773  32.320  33-603 

0.4000  4.96  0.161  1.766  3.001  19-975  31.933  33.103 

0.5000  5.03  0.155  1-593  7-123  13.457  31-923  40.734 

0.6000  5-07  0.143  1.450  6.541  17.573  31.963  42.321 

0.7000  5.09  0.131  1.337  6.149  17-047  32.029  43-303 

0.8000  5.11  0.122  1.254  5-388  16. 716  32.O84  43-937 

0.9000  5-12  0.114  1.196  5-714  16.506  32.127  44.344 

1.0000  5.13  0.109  1.156  5.597  16.369  32.159  44.610 

2.0000  5.15  0.093  1.075  5-374  16.117  32.226  45.110 

Equilibrium  Distribution 

2.3839  5.15  0.093  1.074  5.371  16.113  32.227  45-117 


12 


Table  1-c  Non-equilibrium  charge  distribution  obtained  with 
program  NEDT.   Initial  charge  state  is  6+. 


Ionic  Charge 
Distribution 


X 


ZQX 


Y(l) 


0.0010 

0.0020 

0.0030 

0.0040 

0.0050 

0.0060 

0.0070 

0.0030 

0.0090 

0.0100 

0.0200 

0.0300 

0.0400 

0.0500 

0.0600 

0.0700 

0.0300 

0.0900 

0.1000 

0.2000 

0.3000 
0.4000 
0.5000 
0.6000 
0.7000 

0.8000 
0.9000 
1.0000 
2.0000 


6.00 
5-99 
5.99 
5.9S 
5.96 
5-97 
5.97 
5-96 
5.96 
5.95 
5-91 
5-^7 
5. S3 
5-79 
5.76 

5-73 
5.70 

5.67 
5.65 
5.46 
5.34 
5.2S 
5.23 
5.20 
5.1S 
5.17 
5.16 
5.16 
5.15 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.005 

0.016 

0.032 

0.047 
0.061 
0.072 

0.030 

0.035 

0.039 

0.093 


Non-equilibrium     Initial  Charge  6+ 
Y(2)    Y(3)    Y(4)    1(5)    Y(6) 


0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.001 
0.002 

0.004 

0.006 

0.009 
0.014 
0.019 
0.134 
0.319 
0.510 
0.670 
0.793 

0.331 

0.943 

0.936 
1.016 
1.073 


0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.001 
0.001 
0.006 
0.020 
O.QLL 

0.073 
0.124 

0.130 

0.247 
0.324 
0.409 
1.526 
2.652 

3.513 
4.127 
4.540 
4.313 
5.003 
5.127 
5.210 

5.369 


0.000 
0.003 
0.003 
0.014 
0.023 
O.O33 
0.04.5 
0.053 
0.074 
0.090 
0.332 
0.639 
1.123 
1.626 
2.162 
2.721 
3.291 
3.663 
4.430 
9.163 
12.001 
13.610 
14.543 
15.114 
15. ^67 
15.692 
15.337 
15.932 
16 . Ill 


0.500 

0.933 

1.469 

1.942 

2.407 

2.365 

3. 316 

3. 760 

4.197 

4.627 

3.573 

11.953 

14.839 

17.309 

19.426 

<c  JL .  <c44 

22.306 
24.150 
25.3U 
30.947 
32.235 
32.492 
32.492 
32.436 
32.379 
32.333 
32.300 
32.276 
32.227 


99.500 

99.009 

93.524 

9S.0U. 

97.570 

97.102 

96.639 

96.131 

95.729 

95.232 

91.033 

37.333 

33.939 

30.935 

73.235 

75.349 

73.647 

71.6Z.9 

69.330 

53.221 

52.777 

49.339 

43.116 

47.056 

46.333 

^5.949 

45.665 

'45-477 

45.122 


Equilibrium  Distribution 


2.3143  5.15   0.093   1.074   5.371  16.113  32.227  45.117 
Cross-sections  for  Capture  (IN  10E-16  CM2/MDL) 


N     SC(1)  SC(2)   SC(3)   SC(4)   SC(5) 

1     0.0  1.000   2.000   3.000  4.000 

Cross-sections  for  Loss  (IN  10E-16  CM2/M0L) 

N     SL(1)  SL(2)   SL(3)   SL(4)  SL(5) 

1    11.000  10.000  9.000   3.000  7.000 


SC(6) 
5.000 


SC(7) 


SL(6)   SL(7) 
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ground  state  when  they  undergo  the  next  collision,  and  the 
probabilities  for  charge  exchange  from  and  into  excited 
states  differ  considerably  from  those  for  the  ground  state. 
Then  the  form  of  Eq.  (1)  cannot  be  used.   In  this  work  at- 
tention is  confined  to  very  dilute  gaseous  strippers,  where 
excitations  are  believed  to  have  no  significant  influence. 

An  interesting  quantity  is  the  mean  ionic  charge  I,  the 
center  of  an  equilibrium  distribution.  The  usual  definition 


is 


l 


(2.3) 


where  Y.  are  normalized  charge  fractions  (  £Xj-|  ).  On   the 

other  hand,  highly  asymmetrical  charge  distributions  are  very 

often  observed;  consequently  the  most  probably  charge  i, 

which  is  the  charge  for  which  the  smoothed  distribution  shows 

its  maximum  may  be  defined.   For  symmetrical  distributions, 

i  equals  i,  but  for  asymetrical  distributions,  differences 

i  -  i   0.6  can  be  found.  Further  details  are  discussed  in 
part  V. 

Present  Status  of  Experiments  and  Theories 

A  detailed  review  of  the  investigations  of  oharge  ex- 
change oollisions  prior  to  1965  can  be  found  in  an  artiole 
by  V.  S.  Kikolaev  (1).  Briefly,  the  development  was  as 
follows.  Rutherford  (2),  one  of  the  first  investigators, 
examined  the  passage  of  protons  and    partioles  through 


14 

gases  and  determined  charge  fractions  and  cross  sections. 
Later,  Lassen  (3,  4)  used  fission  fragments  to  study  charge 
exchange  by  heavy  ions.   These  experiments  were  difficult  to 
analyze  because  fission  fragments  show  a  spectrum  of  both 
masses  and  energies.  Moreover,  initial  charge  states  could 
only  be  estimated.  After  the  development  of  particle  ac- 
celerators, intensified  studies  began  with  light  ions  up  to 
argon  inclusive  (see  Reference  1)   In  recent  years  the 
techniques  for  heavy  ion  acceleration  have  been  improved 
and  detailed  experiments  have  been  started  for  heavy  ions 
up  to  uranium  (5-10).   This  is  of  special  interest  in  con- 
nection with  the  design  of  future  heavy  ion  accelerators, 
where  high  ionic  charge  states  are  desired  in  order  to  obtain 
effective  acceleration  to  very  high  energies.   However,  the 
necessity  of  avoiding  collisions  of  accelerated  ions  with 
residual  gas  in  the  beam  line  requires  a  knowledge  of  charge 
exchange  cross  sections.  Although  much  information  is 
already  known,  many  questions  are  still  open  and  more  heavy 
ion  data  must  be  obtained  in  order  to  gain  a  complete  under- 
standing of  stripping  phenomena. 

Theoretical  studies  have  been  done  extensively  by  Bohr 
(11)  (1948),  but  at  that  time  not  much  experimental  data  was 
available  which  could  be  used  to  test  theoretical  results. 
Further  attempts  have  been  made  (12-16),  but  the  situation 
turned  out  to  be  extremely  difficult.  Whereas  the  passage 
of  protons  and    particles  through  H2  and  He  can  be  described 
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very  well,  all  predictions  for  heavier  ions  and  heavier  targets 
at  high  energies  are  very  poor  and  can  hardly  be  used  for 
practical  purposes.   In  all  theories  stationary  models  are 
assumed,  as  well  as  various  approximations,  such  as  the 
Thomas-Fermi  model.   It  seems,  however,  that  an  adequate 
description  can  only  result  from  a  dynamical  model  which  also 
takes  into  account  the  atomic  shell  structure  of  both  the  ions 
and  the  target  atoms.   The  magnitude  of  the  cross  sections 
for  heavy  ions  (KT1*  cm2/atom)  indicates  that  average  impact 
parameters  are  of  the  size  of  the  ions  and  atoms. 

Therefore,  one  should  expect  that  during  an  encounter 
both  the  ion  and  the  target  atom  suffer  considerable  distur- 
bances in  their  electron  shells.   In  addition,  many  shell 
effects  have  been  found  (7,  8)  which  indicate  that  statistical 
models  are  only  of  limited  accuracy.  Finally,  none  of  the 
presently  existing  theories  takes  into  account  the  important 
multiple  electron  loss  cross  sections.   It  is  therefore  ob- 
vious that  improved  models  are  necessary  in  order  to  allow 
reasonable  predictions  for  the  basic  cross  sections  and 
quantities,  such  as  mean  ionic  charges  and  charge  fractions 
under  given  conditions. 
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CHAPTER  III 

METHODS  FOR  NUMERICAL  INTEGRATION  OF  NON-EQUILIBRIUM 
CHARGE  DISTRIBUTIONS 

Whenever  charge  exchange  cross  sections  0^±j   are  to  be 
determined  from  experimental  non- equilibrium  distributions, 
or  where  non-equilibrium  distributions  are  to  be  calculated 
from  given  cross  sections,  numerical  integration  of  the  system 
of  differential  equations  (1)  becomes  necessary.  Although 
various  standard  techniques  exist  for  the  numerical  solution 
of  linear  first  order  differential  equations  with  constant 
coefficients,  it  seems  advantegeous  to  develop  special  pro- 
grams which  satisfy  all  requirements  in  connection  with  heavy 
ion  work.   For  example,  the  determination  of  cross  sections 
from  experimental  charge  distributions  requires  repeated 
evaluation  of  the  system  (1)  with  certain  cross  section  sets; 
thus,  the  integration  routine  should  be  very  fast  with  respect 
to  the  execution  time,  while  the  error  of  the  calculated 
charge  fractions  must  be  within  approximately  1%.      Therefore, 
various  integration  techniques  discussed  below  have  been 
studied. 

It  can  be  assumed  that  initial  conditions  are  given  for 
all  charge  states  i  .  .  .  >  ig: 

Y(i)  =  Y0(i)   for  x  =  xQ  and  i  =  iA.  .  .i£  .        (3.1) 

If  Eq.  (1)  is  used,  ^Yi   can  be  calculated  for  a  small  in- 
crease in  path  length,  ^x: 
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AY(±)   =  ?0\$   Tj(x)  Ax  .  (3.2) 

As  Eq.  (3-2)  can  be  applied  to  all  charge  states  i,  one  can 
find  all  new  components  1%  (  x  +  x) .      If  this  procedure  is 
repeated  the  distribution  Y^(x)  can  be  found  for  all  target 
thicknesses  x  up  to  equilibrium.   Though  it  is  well  known  that 
this  simple  method  is  not  very  accurate,  it  can  be  shown  that 
an  appropriate  choice  of  the  step  width  x  leads  to  errors 
$1   which  are  acceptable  for  our  purpose.  Therefore,  the 
step  width  is  discussed  in  several  aspects. 

Constant  Step  T.'/idth 

All  cross  sections  may  be  arranged  as  a  matrix  M  ~[^±-\l) 
and  the  charge  fractions  at  a  certain  target  thickness  may 
be  arranged  as  a  vector  Y.   Then  for  a  given  YQ  (x  ) , 
1,  (X,+  ax)   can  be  calculated  as  follows: 


\(xQ+   x)  =  (1+  xM)  Y0(xQ)  .  (3.3) 


If  this  procedure  is  repeated  n  times,  one  obtains  the  expres- 
sion 

Yn(x)  =  (1+  xIT)nYo;  x  =  xQ  +  n  x  .  (3.4) 

In  order  to  derive  a  condition  for  Ax  which  guarantees 
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convergence,  the  matrix  M  is  transformed  to  a  diagonaliaed 
matrix  D: 

M  =  SDS"1,  and  D  =  S_1MS  -  [Am J  .  (3-5) 

This  relation  can  be  used  to  rev/rite  3q.  (3-4)  in  the 
following  form: 

Yn  *  S(I+^xD)n  S"^  .  (3.6) 

In  expression  (3-6),  I  +  axD  is  also  a  diagonal  from  and  one 
may  write  the  nonvanishing  elements  as  |-f  &.x\vr\.   It  is 
obvious  that  convergence  is  only  obtained  if  (l+iXAm)   is 
limited.   This  is  equivalent  to  the  condition 

Jl+AxAm|<!     =*AX<|1^FI  (3'7) 

Condition  (3-7)  must  be  fulfilled  for  the  smallest  of  all  the 
negative  eigenvalues  \   .   Unfortunately,  these  eigenvalues 
are  usually  unknown.   But  the  fact  that  the  relation 
causes  one  of  the  eigenvalues  to  assume  the  value  zero  allows 
one  to  estimate: 

^*4W|«|M>|   ,  (3.8! 


where  N  is  the  number  of  charge  states  and  TrD  =■>   Xm  --?(??. 
Sqs.  (3.7)  and  (3.S)  allow  the  estimation  of  Ax,  because  all 
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cross  sections  are  given.   Table  1  shows  three  non-equilibrium 
distributions  which  have  been  calculated  for  three  different 
conditions.  At  x  =  0  ,  all  charge  states  were  assumed  to  be 
zero  except  T,,  Y^,  and  Y6  =  100$  for  the  three  distributions, 
respectively.  The  logerithmic  scale  was  chosen  for  a  better 
demonstration  of  the  range  of  Y  and  x.   It  can  be  seen  that 
a  very  large  number  of  steps  is  necessary  because  x  has  to 
be  chosen  small  compared  with  x  =  1.   On  the  other  hand,  the 
distributions  Y(x)  are  very  smooth  curves  and  change  less 
and  less  if  x  increased,  hence  a  variable  step  width  should 
be  employed. 

Variable  Step  Width 

The  simplest  method  is  to  choose  a  logarithmic  step 
width.  Then,  in  the  logarithmic  plot  (Plate  1)  dx  would 
always  be  represented  by  a  constant  interval.  Thus,  the 
actual  increment  depends  on  x  and  can  be  defined  as 


x  =  x 


f  ,  (3.9) 


where  f  is  a  constant  fraction.  Test  runs  have  shown  that 
values  of  f as  0.03  yield  results  which  are  of  an  accuracy 
better  than  1$.   In  many  cases  the  calculated  fractions  are 
to  be  calculated  in  regular  intervals.  In  Appendix  I,  a 
source  listing  of  a  program  which  integrates  system  (1)  with 
a  logarithmic  step  width  as  described  above  is  given.  A  minor 
modification  of  Ax  was  made  in  order  to  print  the  resulting 
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fractions  ten  times  per  cycle-   In  Appendix  I,  a  flow  chart 
of  the  basic  integration  routine  NEDT  is  also  given.  In  this 
program,  the  number  of  integration  steps  per  cycle  can  be 
controlled  by  the  input  data.  The  required  input  information 
is  also  listed  in  Appendix  I. 

Accuracy  and  Execution  Time 

Various  test  runs  have  been  made  in  order  to  find  a 
suitable  constant  f  as  defined  by  Eq.  (.3.9).   In  general,  if 
f  is  decreased,  the  accuracy  will  become  better  at  the  cost 
of  execution  time.  On  the  other  hand,  for  valxies  of  f  smaller 
than  a  certain  minimum,  computational  errors  will  not  allow 
any  improvement.   In  Plate  2,  two  charge  states  have  been 
selected  where  the  errors  AY  have  a  maximum.  Due  to  the 
fact  that  both  at  small  and  large  values  of  x  the  uncertain- 
ties become  negligible  (see  part  "Equilibrium"  below),  these 
maximum  errors  occur  somewhere  between.   Plate  2  indicates 
that  300  steps/cycle  gives  an  accuracy  better  than  about  1%. 
This  is  acceptable  in  many  cases  where  measured  charge  frac- 
tions with  uncertainties  of  5~10?o   are  to  be  compared  with 
calculated  fractions.  The  exact  values  in  Plate  2  have  been 
determined  by  means  of  the  Runge-Kutta  method  (see  below). 

Equilibrium 

It  can  be  shown,  that  the  equilibrium  distribution  is 
unique  for  given  cross  sections.   If  the  index  n  in  Eq.  (3.4) 
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becomes  infinitely  large,  which  corresponds  to  equilibrium, 
obtains 

Yoo  =  ( 1+  xM)  Y«,  *  MY^O  .  ( 3  •  10 ) 

If  this  result  is  compared  with  the  differential  Eq.  (1), 
which  can  be  written  as  YT  =  MY,  it  is  apparent  that  Y'  =  0 
is  indeed  the  condition  for  equilibrium.  Thus,  whereever 
this  integration  method  yields  equilibrium  it  must  be  the 
correct  one  although  charge  fractions  at  smaller  x  may  have 
large  errors.   It  should  be  noted,  however,  that  due  to 
computational  errors  the  condition  Y'  =  0  can  not  be  fulfilled 
exactly.   If  the  integration  procedure  is  not  interrupted, 
the  fractions  Yi  will  not  remain  constant  but  will  continue 
to  change,  and  may  assume  unphysical  values.   In  this  program, 
a  test  for  equilibrium  is  made.   If  the  total  change  of  all 
fractions  becomes  smaller  than  a  certain  minimum  amount  the 
integration  procedure  is  stopped. 

Runge-Kutta  Integration 

Other  methods,  especially  the  Runge-Iutta  integration 
techniques,  have  also  been  studied.  Suitable  application  of 
the  well  known  Runge-Kutta  method  gives  extremely  accurate 
results.  In  fact  one  could  not  produce  fractions  with  obser- 
vable errors  without  losing  convergence.   Therefore,  the 
results  obtained  with  this  method  could  be  used  for  a  comparison 
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with  the  simple  integration  techniques  described  above.  The 
only  disadvantage  is  that  the  Runge-Kutta  method  requires 
much  more  execution  time.   Test  runs  indicated  a  time  ratio 
of  at  least  1:4. 
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CHAPTER  IV 
DIRECT  CALCULATION  OF  EQUILIBRIUM  CHARGE  DISTRIBUTION 

It  has  been  shown  in  Chapter  III  that  a  unique  equili- 
brium distribution  can  be  derived  from  a  given  set  of  cross 
sections.  Although  the  integration  programs  described  in 
Chapter  III  also  give  equilibrium,  a  more  direct  method  can 
be  of  advantage  for  various  reasons.  Application  of  the 
integration  technique  requires  the  assumption  of  initial 
conditions  for  the  charge  states,  which  do  not  effect  the 
equilibrium.   In  addition,  the  execution  time  is  unneces- 
sarily long  if  non- equilibrium  fractions  are  calculated  but 
not  desired.  The  most  direct  way  for  a  calculation  of  the 
equilibrium  charge  distribution  is  to  solve  the  equation 

MY^=  0  .  (4.1) 

In  fact,  the  whole  calculation  then  is  reduced  to  solving 
a  linear  inhomogenous  system  of  equations.  Equation  (4.1) 
turns  out  to  be  an  inhomogeneous  system  because  the  relation 
2Yoa=  1  kas  to  ^e  used  to  reduce  (4.1).   Otherwise  only 
trivial  solutions  would  exist.   If  one  eliminates,  for  ex- 
ample, YkQaby  Yk  =1-  £  Yj_  00  ,  one  gets  the  form 

iVk 

1 

MYoa=  C  ,  (4.2) 

with  the  elements 


26 


and  Ci  =  -Mik  -  ^  . 

The  remaining  problem  is  equivalent  to  finding  the  inverse 
matrix  M~  .   Then  the  desired  fractions  are  obtained  by  the 
products 


Y  =  M"1  C  .  (4.3) 


The  common  standard  technique  is  known  as  row  operation.   In 
this  method,  the  diagonal  elements  of  M  are  successively- 
changed  to  unity,  whereas  the  off -diagonal  elements  become 
zero.  Combininb  M  and  C  in  a  partition  matrix,  the  operations 
can  be  symbolized  as  follows: 

[M,  C]  =>    [i,  C]  ,  (4-4) 

where  the  elements  of  Cf  become  the  desired  fractions  Y 
The  remaining  value  YK«,  is  then  obtained  f rom  X^k  =l~Z  Yi»  • 
This  technique  is  very  effective  and  takes  a  minimum  of  com- 
putational time.  As  all  diagonal  elements  of  M  are  always 
different  from  zero,  no  complication  can  arise  and  no  measures 
have  to  be  taken  in  order  to  correct  for  an  improper  matrix. 

The  detailed  source  listing  for  the  program  EQUDIS  is 
also  given  in  Appendix  I  as  a  subroutine,  including  a  brief 
flow  chart.  As  can  be  seen  from  Table  1,  there  is  perfect 
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agreement  between  the  two  equilibrium  distributions  calculated 
from  the  integration  program  NEDT  and  the  routine  EQUDIS. 
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CHAPTER  V 
ANALYSIS  OF  EQUILIBRIUM  OF  THE  BR  AND  I  ION  METHOD 

In  this  part,  some  of  the  equilibrium  charge  distribu- 
tions are  analyzed  which  have  been  measured  at  Oak  Ridge  (7). 
Our  interest  is  not  only  in  the  mean  ionic  charge  i  but  also 
in  the  shpae  of  the  distribution;  semi-empirical  formulas  for 
i  can  be  tested  as  well  as  possible  relations  between  asym- 
metries of  the  equilibrium  distributions  and  multiple 
electron  loss  processes. 

The  distributions  have  been  measured  for  Br  ions  at  20 
MeV  and  for  I  ions  at  15  and  60  MeV  in  dilute  gases  of  Hv?, 
He,  Np,  Ne,  Ar  and  Kr.  A  common  test  for  a  Gaussian  dis- 
tribution Y(i)  is  to  plot  the  logarithim  of  the  ratios 
T(i)/T(i+1)  versus  the  charge  i.  A  Gaussian  shape  will  result 
in  a  straight  line  and  deviations  can  be  seen  very  easily 
(7).   The  advantage  of  this  technique  is  that  an  incorrect 
normalization  of  the  fractions,  perhaps  due  to  missing  charge 
states,  cannot  cause  errors  since  only  ratios  are  decisive. 
On  the  other  hand,  all  distributions  are  clearly  asymmetrical 
and  one  may  try  to  find  a  better  approximation  for  Y(i)  than 
is  given  by  a  Gaussian  distribution.   Therefore  one  tries  to 
fit  the  equilibrium  charge  distrib\ition  to  a  generalized 
Gaussian, 

T1  =  Y  exp  (-(~)n)  ,  (5.1) 
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where  Y  is  the  maximum  intensity  of  a  smoothed  distribution 
which  occurs  at  the  most  probably  charge  i.  Of  course,  1 
need  not  necessarily  be  an  integer.  The  assumption  (5.1) 
is  arbitrary  and  other  functions  could  be  used.   But  (5-1) 
implies  a  Gaussian  for  n  =  2  and  one  can  expect  to  obtain 
characteratic  relations  between  the  exponent  n  and  the  width 
parameter  J"  for  the  various  fitted  distributions. 

Before  one  can  find  the  best  values  of  n  and  f   for  the 
distributions  in  question  by  means  of  a  least  square  tech- 
nique, Y  must  be  determined.  As  the  fit  does  not  depent  too 
strongly  on  Y,  an  approximation  for  Y  is  sufficient.  There- 
fore, one  defines  Y  as  the  maximum  of  a  parabola,  which  is 
determined  by  the  three  largest  charge  fractions: 


?.  Y.=|  (i  _  I)2    .      '  (5.2) 

If  the  coordinates  i±,"I±)   of  the  three  largest  charge  frac- 
tions  are  inserted  in  Eq.  (5-2),  the  three  unknowns  Y.  P  and 

A 

i  are  obtained: 

a       P  ,     «*  2 
Y  =  Y1+  2(^-1)   , 

P  -  2Y2  -  Y1  -  Y3  ,  (5.3) 

l  =  Y2~Y1  +'i,  +  i  . 


In  order  to  apply  the  least  square  methods  for  a  straight 
line,  Eq.  (5.1)  must  be  transformed  into  a  linear  equation: 


Y  =  nx  +  a  , 


(5-4) 
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where   Y«  In (In  (Y/Y,  ) )  , 
x  -  In  | Ai(  , 

a  =  -nlntf  . 

Applying  the  standard  formulas  for  the  best  fit  of  a  straight 
line,  one  obtains  for  the  best  values  n  and  a: 

n   -  (Ixf-lY^Ix,-  IxiY^/D  ; 

I-  (Nlx^-lx.Y^/D  ;  ^^ 

D  =  N£xf-(  Ix^2  , 

where  all  N  experimental  points  are  to  be  taken  into  account 
in  the  summation.   The  most  probable  error  of  n  and  a  can  also 
be  derived  as 

B  -  (5Lx2/D)Vo  , 

a-  (b/D)*Po  ,         '  {5.6) 

where  N  is  the  number  of  experimental  points,  and  f0  is  given 
by 

fo  ■=  U(y±-Z-1Zxi)2/(iL-2)']  *  (5-7) 

Application  of  this  method  implies,  that  the  two  wings  of  the 
equilibrium  distributions  are  examined  separately.   This  pro- 
cedure seems  to  be  meaningful  because  of  the  asymmetries  which 
will  result  in  different  values  of  n  and  ;   for  the  left  and 
the  right  sides. 
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Results  and  Discussion 

The  results  of  the  various  fits  are  shown  in  Plates 
3;  to  13  and  in  Table  2.  Two  different  fits  have  always  been 
attempted.  First,  n  and  J  have  been  determined  according  to 
the  method  described  above.   Secondly,  an  independent  deter- 
mination of  the  width  fa.  was  made  under  the  assumption  of  a 
Gaussian  with  n  -  2. 

V/henever  one  side  contained  more  than  three  measured 
points,  the  fit  was  made.  In  most  cases,  where  values  of  n 
and  r  could  be  determined  for  the  left  side,  the  result  of 
H  was  close  to  2.  This  means,  that  the  left  sides  of  all 
equilibrium  distributions  are  mostly  Gaussians,  independent 
of.' the  ion  and  the  stripper.  Only  for  I  at  60  MeV,  somewhat 
larger  values  of  n  were  found,  indicating  that  the  decrease 
Y.  to  smaller  charges  i  is  steeper  than  for  a  Gaussian. 

In  all  cases,  the  analysis  of  the  right  sides  of  the 
distributions  resulted  in  values  n  smaller  than  2,  which  in- 
dicates a  less  rapid  decrease  than  Gaussian.   In  Br,  n  ■  1.3 
was  found  for  all  targets  whereas  in  iodine  IS  shows  some 
dependence  on  the  stripper  (see  Table  2).  The  width  para- 
meter for  the  right  sides  is  smallest  for  Br,  increases  for 
r  at  15  MeV,  and  is  largest  for  I  at  60  MeY.  The  width  does 
not  depend  too  much  on  the  target  except  in  the  case  of  H2, 
where  smaller  values  of  f    are  always  found.   The  results 
are  given  only  for  comparison.  As  the  corresponding  fit, 
especially  for  the  right  sides,  is  very  bad  under  the 
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Table  2  Distribution  parameters  and  mean  ionic  charges  for 
Br  and  I  ions  in  various  gases  as  obtained  from 
the  fit  procedures. 


Left  Side      Right  Side 
n  n 


Br  20  MeV  H2  2.33  1.66  1.73  I.96  1.82  1.81  6.17  6.20 

He  2.01  1.74  1.74  1.31  1.81  2.12  6.88  7.11 

N2  2.01  1.68  1.70  1.26  1.99  2.30  7-02  7.37 

Ne  2.00  1.79  1.80  1.28  2.00  2.19  6.87  7-06 

A  2.00  1.73  1.71  1.33  1.93  2.00  6.92  7.14 

Kr  2.01  I.64  1.62  1.30  1.95  2.15  7-01  7.16 

I  15  MeV  H2  2.16  1.88  1.90  1.67  2.10  2.21  4-44  4-5S 

He  —  —  1.36  2.46  2.65  4-93  5-54 

N2  —  —  —  1.46  2.60  2.69  5.00  5-69 

Ne  —  —  —  1.40  2.61  2.68  5.O4  5.76 

A  —  —  1.32  2.45  2.64  4-96  5.64 

Kr  —  —  —  1.24  2.34  2.68  4.77  5.45 

I  60  MeV  He  2.06  2.14  2.18  1.682  3.20  2.72  11.15  11.58 

N2  2.30  2.46  2.50  1.82  3.20  3.20  12.5S  13.09 

Me  2.17  2.45  2.65  1.5S  3-10  3-32  11.54  12.07 

A  2.37  2.50  2.54  1.86  3.24  3.23  12.99  13.48 

Kr  2.35  2.50  2.53  1.86  3-30  3.30  12.77  13-13 
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assumption  n  =  2,    ^should  not  be  weighted  too  much. 
Plates  3  to  13  give  a  better  idea  about  the  quality  of  the 
various  fits. 

As  only  a  few  measurements  have  been  done  for  heavy  ions 
in  that  energy  range,  it  is  difficult  to  make  a  comparison 
with  other  results.   In  gaseous  targets  the  trend  seems  to 
be  as  follows:   for  light  ions,  all  measured  charge  distri- 
butions show  a  nearly  perfect  Gaussian  shape  (1)  in  all 
stripper  gases.  The  heavier  the  ion,  the  greater  the  devia- 
tion from  a  Gaussian;  in  all  cases  the  higher-  charge  states 
have  more  intensity  than  compared  with  a  Gaussian.   The 
stripper  apparently  has  some  influence:   in  H2  and  He,  the 
distributions  are  more  symmetrical  and  in  very  heavy  strip- 
pers the  largest  asymmetries  can  be  found.  This  is  also 
obvious  from  the  differences  l  -  i,  as  listed  in  Table  2. 

One  could  try  to  ascribe  this  asymmetry  to  the  occurrence 
of  multiple  electron  loss.  In  fact,  if  one  increases  multiple 
electron  loss  probabilities,  the  equilibrium  fractions  tend 
to  shift  to  higher  charge  states.   This  is  consistent  with  the 
result  that  multiple  electron  capture  is  very  unlikely,  and 
therefore  no  corresponding  trend  to  smaller  charge  states 
exists.  On  the  other  hand,  even  for  light  ions  such  as  argon, 
where  Gaussian  distributions  have  been  measured,  multiple 
electron  loss  can  be  significant.   Obviously  the  equilibrium 
distribution  is  a  result  of  the  competing  processes  for 
electron  capture  and  loss,  C~±y         But  these  cross  sections 
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depend  on  the  electronic  shell  structure  of  both  the  ion  and 
the  target.  Therefore  one  must  expect  different  shapes  for 
the  distributions,  depending  on  the  ion,  their  energy,  and 
the  target.  Also,  shell  effects  will  occur  (7).   It  is 
somewhat  surprising  that,  nevertheless,  very  similar  equilib- 
rium distributions  can  be  observed;  for  example,  for  one  ion 
in  one  target  in  a  very  wide  energy  range.  Although  the 
cross  sections  which  determine  the  distribution  depend  strongly 
on  the  energy,  it  seems  as  if  the  distribution  did  not  change 
its  shape  but  is  simply  shifted  to  higher  charge  states  as 
the  energy  increases  (7). 

We  conclude  that  empirical  formulas  for  the  equilibrium 
distributions  of  fast  ions  can  be  found,  but  due  to  the 
specific  atomic  structure  of  both  ions  and  target,  no  simple 
relation  Y(i)  can  be  given  which'  holds  in  a  general  sense. 
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Appendix  I. 

Source  listing  of  programs  NEDT  and  EQUDIS,  controlled  by- 
main  routine,  a  list  of  input  parameters  and  a  brief  flow 
chart . 

MAIN  PROGRAM 

COMMON  CS  (30,30),  JA(30),JE(30),ZF(9),IO,IA,IE,NC,NL,IPO, 

YP,X?,X,NP 

1  FORMAT (1H1,A6,4S,2CHNUC LEAR-CHARGE  2  -  , 14 , 11X , 17KAR0MIC 
IM1ASS  M  =,H,//,11X,  20HKIN.  ENERGY  E  =,F6.2,26H  MEV 
2VEL0CITY  V/VO  =,F7.3,//,9H  STRIPPER,A8,14H  -GAS  10  =,13, 
31H+,11X,13HCHAR  GE  STATES, 13, 4H+  TO,13,lH+,//,31H  MULTIPLE 
4STRIPPING  NC  =,I3,12X,4HNL  =,12,///} 

2  FORMAT (1H-,  35HCROSS-SECTION  FOR  CAPTURE  (IN  10E-,12, 
19H  CM2/M0L)) 

3  FORMAT  (1H0,6X,13(6X,I2,1X)) 

L       FORMAT (1H+,6H  N  ,13(3X,6HSC(   ) ) ,/) 

5  FORMAT   (10(F6.3,2X)) 

6  FORMAT  (1H  ,13,3X,13F9-3) 

7  FORMAT (1K-,35KCR0SS-SECTI0NS  FOR  LOSS  (IN  loE-,12,9K 
1CM2/M0L)) 

8  FORMAT (1H+,6N  N  ,13(3X,  6HSL(   )),/} 

9  FORMAT  (A4,12,13,F6.2,F6.3,A4,813,3F7.5) 
10   DATA  IRE,IWR/l,3/ 

DATA  IST,IG0/4HST0P,4HG0  / 

14  READ  (IRE, 9)  IN,IZ,IB,E,V,STR,IA,IE,N?,IPO,NC,NL,LO,IQ, 
1X,YP.XP 

IF   (IPO.EQ.O)  IP0=l6 
IF  (IN.EQ.IST)  STOP 
IF  (IN.EQ.IGO)  GOTO  98 

15  WRITE  (IT,7R,1)  INjIZjIB.E.VjSTRjIO.IAjIE^^jNL 
WRITE(IWR,2)  IPO 

WRITE (  r.7R,3)  (1, 1=1  A,  IE) 

WRITE(IV/R,4) 

IF  ((IE-IA).GT.L#)   STOP 

DO  16  I=IA,IE 

DO  16  J=IA,IE 

16  CS(I,J)=0. 
DO  18  M=L,MC 
K-N+LA 
I1=IE 

IF  ((IE-IAJ.GT.9)  Il=IA+9 

READ  (IRE,  5)  (ZF(L),L=1,N),(CS(I,I-U),I=K,I1) 
IF  ((IE-IA).LE.9)  GOTO  18 
12=11+1 

READ   (IRE.  5)  (CS(I,I-N),I=I<2,IE) 
18   WRITE(IWR,6)  N,(ZF(L),L=1,N),(CS(I,I-N),I=K,IE) 
WRITE(r.7R,7)  IPO 
WRITE(IWR,3)  (I,I=IA,IE) 
WRITE (IWR,g) 
■  DO  20  N=1,NL 
K-IE-H 
I1=K 
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IF((K-IA).GT.9)  I1-IA+9 

READ (IRE, 5)  (CS(I,I+N),I-IA.I1) 
IF  ((K=IA).LE.9)  GOTO  20 
12=11+1 

READ (IRE, 5)  (CS(I,I+N),I-I2,K) 
20   WRITE  (IWR,6J  N,  (CS(I,I+N),I-IA.K) 

DO  25  I=IA,IE 

KA=MACO(I-NC,IA) 

IE=MINO(I+NL,IE) 

SUI  1=0 . 

DO  24  K=KA,KE 

24    sum=sum-cs(i,k) 

CS(I,I)-SUM 

JA(I)=MAX0(I-NL,IA) 
25  JE(I)=:HNO(I+NC,IE) 
93   IF  (IO.NE.O)  CALL  NEDT 

IF  (IQ.NE.O)  CALL  EOUDIS 

GOTO  14 
99   STOP 

END 

SUBROUTINE  NEDT 

COMMON  CS(30,30),JA(30),JE(30),ZF(9),IO,IA,IE,NE,NL,IPO, 
YP,XP,X,NP  DIMENSION  Y(10),YS(30)   »>»»>> 

1  FORMAT (1H1, 43 KICNIC  CHARGE  DISTRIBUTION  NON-EQUILIBRIUM, 
112X,  14HINITIAL  CHARGE, 13, 1H+,9X.9HX  IN  I0E+,12,26H 
2M0L/CM2  Y  IN  PERCENT, 17, 6H  STEPS) 

2  FORMAT  (IH-,15X,13(6X,12,i:--)) 

3  FORMAT (IH+,4X,1HX,6X,4H  ZQX,13 (4X,5HY(   )),/) 

4  FORMAT  (IK  ,FS.2I.,F7.2,2PF9.3,12F9!3) 

5  FORMAT (1H0,  26HEQUILIBRIUM  DISTRIBUTION   ,/) 

6  FORMAT(5X,FS.3,6X,6(3X,F7.3)) 

7  FORMAT(3(3X,?7.3)i 

8  FORMAT ( IK- ,  45HEQUILIBRIUM  DISTRIBUTION  COULD  NOT  BE  REACHED) 
DATA  IWR,IPU/3,2/  ; 
IF  (NP.EQ.O)  NP=10 

PN=1. /FLOAT (NP) 
IF  (NP.EQ.10)  PN-.l 
WRITE ( IWR , 1 )  10 , I n0 , NP 
WRITE ( IWR, 2)  (I,I=IA,IE) 
WRITE (IWR, 3 ) 
IF  (X.EQ.O.)  X=.001 
DO  10  I=IA,IE 
10    Y(I)=X*CS(IC,I) 
Y(I0)=Y(IC)+1. 
ZQX=IO 
YM2=99 . 
DO  16  NN-1,6 
DX»PN*X 
DO  16  11=1,9 

WRITE(IWR,4)  X,ZQX,(Y(I),I-IA,IE) 

IF  (YP.EQ.C.)  GOTO  12 
IF  (X.LT.XP)  GOTO  12 
DO  11  I=IA,IE 
YS(I)0. 
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IF   (Y(I).LT.YP)   GOTO  11 
YS(I)=Y(I)*100. 

11  CONTINUE 
Il=IA+5 

WRITE(IPU,6)  X,(YS(I),I=IA,I1) 
IF  (IS.LE.I1)  GOTO  12 

12+11+1 

WRITE (IPU,7)  (YS (I) ,1=12, IE) 

12  DO  16  M=1,NP 
X=X+DX 
ZQX=0 . 
YI-n>YL!2 
IM2-0 . 

DO  H  I+IA,IE 

YS(I)=0. 

JAA=JA(I) 

JEE=JE(I) 

DO  13  J=JAA,JEE 

13  YS(I)=YS(I)+CS(J,I)*Y(J) 

14  YM2=YM2+ABS  ( YS  ( I J ) 
DO  15  I+IA,IE 
Y(I)«=Y(I)+YS(I)*DX 

15  ZQX=ZQX+ Y ( I ) *FLOAT ( I ) 
IF  (BI2.GY.YIH)  GOTO  17 

16  CONTINUE 

17  IF  (YM2*DX.LE..01)  GOTO  18 
WRITE (IWR,g) 

RETURN 

18  WRITE(IWR,5) 

WRITE(IWR,4)  X,ZQX,(Y(I),I=IA,IE) 

RETURN 

END 

SUBROUTINE  EOUDIS 

DOIOT  DS(30,3O),JA(3O),JE(3O),ZF(9),IO,IA,IE,NC,NL,IPO, 

YP  j XP  ,  X , NP 

DIMENSION  A (30, 30) 

1  FORMAT  (lH-,rlHEXACT  EQUILIBRIUM  DISTRIBUTION  AS  FOLLOWS,/) 

2  F0RMAT(1H-,32HDIAGCNAL-ELEMENT  IS  ZERO  STOP,/) 

3  FORMAT  (13X,3HZQX,13(6X,I2,1H+)) 

4  FORMAT  (1H   ) 

5  FORMAT  (9X,F7.2,PF9.3,12F9.3) 
DATA  r.7R/3/ 

IEE=IE-1 
DO  11  I=IA,IEE 
A(I,IE)=CS(IE,I) 
DO  11  J-IA,IEE 
11   A(I,J)=CS(J,I)-CS(IE,I) 
WRITE  (IWR,1) 
DO  15  I+IA.IEE 
TBMF-A(I,lj 
IF(TEMP.NE.O)  GOTO  12 
WRITE  (IUR,2) 
STOP 
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12  DO  13  J-IA,IE 
IF(J.LT.I)  GOTO  13 
A(I,J)-A(I,J)/TEMP 

13  CONTINUE 

DO  15  KKEA.IEE 

TEMP=A(K,l} 

IF (K. SO. I)  GOTO  15 

DO  14  J-IA,IE 

IF(J.LT.I)  GOTO  14 

A(K,J)=A(K,J)-TEMP*A(I,J) 

14  CONTINUE 

15  CONTINUE 
ZQ=0. 
A(IE.IE)-1. 

DO  16  I=IA,IEE 
A(I,IE)=-A(I,IE) 
ZQ=ZQ+A(  I,IE*FLOAT(  I) 

16  A(IE,IE)=A(IE,IE)-A(I,IE) 
ZQ=ZQ+A  ( IE ,  IE*FLOAT  ( IE ) 
WRITE  (B7R,3)  (I, MA, IE) 
WRITE  (IWR,4) 

WRITE  (IWR,5)  ZQ,(A(I,IE),I  IA,IE) 

RETURN 

END 
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INPUT  PARAMETERS 


CONTROL  CARD:   (format  9) 


IN 

IZ 

IB 

E 

V 

STR 

IA 

IE 

NP 

IPO 

NC 

NL 

10 

IQ 


X 


YP 


TV 


ion,# 

nuclear  charge 
mass 

energy  in  MeV 
velocity  in  VQ 
stripper 

smallest  charge  state   0 
highest  charge  state 

NP  x  50  steps  per  cycle  (if=0,  NP  =  2) 
unit  for  CS  and  X,  (if  -  0,   IPO  -  16) 
maximum  value  of  multiple  capture  (  =  0) 
maximum  value  of  multiple  loss  (=  0) 
initial  charge  (if  =  0,  no  CALL  NEDT) 
=  0.  =  exact  equilibrium  charge  dist.  will  be 
calculated. 

=  0.  «=  no  extra  calculation  equilibrium  charge 
distribution, 
starting  target  thickness  (if  0,  x  -  0.001  x 
10-IP0  atom/cm2) 
=  0  =  Y  are  punched 
=0    Y  are  not  punched 
if  cards  punched,  then  only  if  X   XP 
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CONDITIONS: 

1.  IA   0 

2.  IE-IA   12 

3.  IA  10  IE 

#IN  =  STOP   program  stop 
IN  =  GO   CALL  NEDT  for  another  initial  charge  state  and 
same  CS  set 
*not  used  for  calculation 


START 
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SUB.NEDT 


READ  DATA 
ARRANGE  M 


SUB.EQUDIS 


DEFINE  Y, 


*  INCREAMENT  A  X 


NO 


CHANGE  AY 


NEW  FRACTIONS  Y 


CHECK  STEP  NUIISER 


PRINT  Y 


YES 


REDUCE  H  TO  M 


M  ILL  CONDITIONED  1 


NO 


DIAGONALIZE  M 
BY  ROW  OPERATION 


CALCULATE  MEAN 
IONIC  CHARGE  I 


PRINT 
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The  processes  of  electron  capture  and  loss  of  fast 
heavy  ions  passing  through  matter  are  described  and  the 
status  of  relevant  experimental  and  theoretical  work  is 
reviewed.  Programs  for  the  numerical  calculation  of  ionic 
charge  states  in  no n- equilibrium  and  equilibrium  are  pre- 
sented, which  can  be  used  for  given  sets  of  charge  exchange 
cross  sections.   Experimental  equilibrium  distributions  of 
Br  and  I  ions  at  15,  20  and  60  MeV,  stripped  in  various 
gases,  are  analyzed  with  respect  to  their  form. 


